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1 .  INTRODUCTION 


The  purpose  of  this  paper  is  to  introduce  an  efficient  meth¬ 
od,  which  is  quite  general  and  easy  to  use,  to  the  solution  of  the  gun 
dynamic  problems,  to  describe  some  of  the  features  associated  with  a 
finite  element  computer  program,  and  to  present  some  initial  results. 


The  basic  concept  of  unconstrained,  adjoint  variational  for¬ 
mulation  for  linear  problems  was  described  in  an  earlier  paper  [1]. 

Its  advantage  over  constrained  methods  in  obtaining  approximate  solu¬ 
tions  has  been  demonstrated  for  both  conservative  (self-ad joint)  and 
unconservative  (nonself-adjoint)  problems  [2] .  In  comparison  with 
Galerkin  procedure,  the  unconstrained,  adjoint  variational  formulation 
has  a  further  advantage  in  the  freedom  of  selecting  shape  functions 
which  have  less  requirement  on  differentiability  and  which  are  not 
required  to  satisfy  any  of  the  end  conditions.  The  same  concept  was 
extended  to  solution  formulation  of  initial  value  problems  [3] .  In 
view  of  the  generality  of  this  approach  and  its  easy  adaptability  to 
finite  element  discretizations,  it  appears  to  be  quite  attractive  in 
seeking  solutions  to  the  complicated  problems  associated  with  the 
dynamics  of  gun  systems. yv\ 

The  basis  of  the  present  formulation  for  more  general  cases 
has  been  given  in  a  previous  paper  [4] .  The  special  problem  of  a  uni¬ 
form  gun  tube  is  treated  here  for  demonstrated  purposes.  The  differ¬ 
ential  equation  and  initial  boundary  conditions  are  given  in  Section  2. 
An  unconstrained  variational  statement  which  is  equivalent  to  the  given 
governing  equations  is  stated  in  Section  3.  The  details  of  some  of  the 
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special  features  on  finite  element  discretization  are  described  in 
Section  4.  Finally  some  preliminary  results  are  presented  in  Section 

5. 

2.  GOVERNING  EQUATIONS 

The  motion  of  a  gun  tube  modeled  by  the  lateral  deflection 
of  an  Euler-Bemoulli  beam  is  shown  in  Fig.  1.  The  differential  equa¬ 
tion  in  nondimensional  form  is 

y""  +  (_p  +  g  sin  a)[(l-x)y']'  ♦  y2y 

=  -Py"H(x-x)  -  y2ro[32t2y"  ♦  2St£'  +  y]6(x-x) 

- (gm  cos  o)6(x-x)-g  cos  a  (1) 


where 

P  =  iTR2p 

y  =  y(x,t),  the  tube  deflection 

x  =  spatial  axis  along  the  tube's  length,  0  £  x  £  1 
t  =  time  axis,  0  £  t  <_  T,  T  is  the  time  limit  of  interest 
a  =  elevation  angle 
m  =  projectile  mass 

3  =  acceleration  of  the  projectile,  assumed  to  be  constant 
p  -  bore  pressure,  assumed  to  be  constant 
g  =  gravitational  acceleration 
x  =  (1/2) St2,  projectile  position 
H(x)  =  the  Heavenside  step  function 
6(x)  =  the  Dirac  delta  function 
R  =  tube  inner  radius 
Y  =  c/T 

The  constant  c  is  defined  by  (pA£.V(EI)) where  p,  E  are  density. 
Young's  modulus  of  the  tube  material;  i.  A,  I,  the  length,  area,  and 
second  moment  of  cross-section  of  tube,  respectively  (see  Ref.  [5]). 
Also  in  Eq.  (1),  a  prime  (')  denotes  a  differentiation  with  respect  to 
x  and  dot  (•)*  a  differentiation  with  respect  to  t.  The  derivation  of 
this  equation  and  the  end  conditions  which  follow  have  been  previously 
given  [4,5]  and  will  not  be  repeated  here. 

The  initial  condition,  or  more  appropriately,  the  end  condi¬ 
tions  in  time  are 

y(x,0)  =  0  ;  y(x,l)[l  ♦  m6(j  S-x)]  ♦  k7[y(x,0)  -  Y(x)]  «  0  (2) 
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Fig.  1  A  schematic  drawing  of  the  problem  configuration 


where  the  constant  ky  is  introduced  in  conjunction  with  the  uncon¬ 
strained  variational  formulation  (Section  3)  so  that  if  one  takes  ky 
to  be  infinite,  the  initial  displacement  y(x, 0)  is  forced  to  be  iden¬ 
tical  to  the  prescribed  shape  Y(x).  Based  on  similar  reasonings  the 
boundary  conditions  have  been  shown  to  be  the  following: 

y"(0,t)  -  k2y'(0,t)  -  0  (3a) 

y"(l,t)  ♦  k4y'(l,t)  -  0  (3b) 

y"'  (0,t)  ♦  kiy(0,t)  ♦  (-P  ♦  g  cos  )y'(0,t) 

♦  Py»(0,t)H(j  8t2)  ♦  m6V(0,t)5(y  ^  =  0  (3c) 

and 

y"'  (l,t)-kj(l,t)*Py’  (l.t)H(^  Ut2-l)*meV  (l,t)6(y  6t2-l)  -  0  (3d) 

where  k^,  i  «  1,  2,  3,  4,  are  the  appropriate  elastic  spring  constants 
at  the  supports. 
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3.  UNCONSTRAINED  VARIATIONAL  STATEMENT 

Through  integrations-by-parts,  it  is  straight  forward  to 
show  that  the  following  variational  statement  is  equivalent  to  the 
differential  equation  and  end  conditions  stated  in  Section  2. 


12 


51  s  C«DV  =  I  («I,)V  -  1  (WJ  -  0 
y  i=i  1  y  j*i  J 

with  1  1  _  11 

I  »f  /  y"y*"dxdt  ;  I“(P-g  sin  a)/  /  (l-x)y'y*'dxdt 

1  0  0  1  0  0 


(4) 


11..  .11 

I  =-y2 3J  /  yy*dxdt  ;  I  «-pJ  /  y'y*’H(x-x)dxdt 


0  0 


0  0 


.  1  1  11 
I  *-pJ  J  y’y*'fi(x-x)dxdt  ;  I  =-mS2y2J  /  t2y'y*'5(x-x)dxdt 


0  0 


0  0 


11  -  .  1  1 

I  =-mg2y2J  /  ty,y*,6' (x-x)dxdt  ;  I  =2m8y2/  /  ty'y*$(x-x)dxdt  (5) 

7  O  A  a 


0  0 


0  0 


1  1"  1  1. 

I  =-my2/  /  yy*6(x-x)dxdt  ;  I  =-my2/  /  yy*6(x-x)dxdt 

9  0  0  10  0  0 


1irJL{kiy(0,t)y*(0,t)  + 


+  kj(l.t)y*(l.t)  +  k  y*  (l,t)y*'(l,t)dt 
3  4 


and 


I  -k  /  y(x,0)y*(x,l)dx 
12  7  o 

1  1 

J  =  -g  cos  a /  /  y*dxdt 

1  0  0 
1  1  .  - 

J.  »  -gm  cos  a /  /  y*6(x-x)dxdt 

1  1  0  0 

J-  -  k _/  Y(x)y*(x,l)dx 

3  7  o 


(6) 
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The  variational  statement  of  Eq.  (4)  will  serve  as  the  basis 
of  our  finite  element  solutions. 

4.  FINITE  ELEMENT  DISCRETIZATION 

The  purpose  of  the  discretization  is  to  enable  one  to  write 
the  variational  statement  of  Eq.  (4),  which  is  a  functional  of  contin¬ 
uous  functions  y  and  y*,  etc.,  in  the  form  of  a  matrix  equation 

6Y*T  K  Y  =  5Y*T  F  (7) 

where  X,  X*  are  the  "global"  generalized  coordinates  vectors.  &  is  the 
"stiffness"  matrix,  and  £  the  "force"  vector.  These  terminology  are 
patented  after  the  static  structural  analysis,  but  they  do  not  neces¬ 
sarily  have  the  physical  meanings  of  those  adjectives  attached  to  them. 
Since  the  variational  statement  associated  with  Eq.  (7)  is  unconstrain¬ 
ed,  the  equation  leads  directly  to 

K  Y  =  F  (8) 

*y  *«• 

which  can  be  solved  for  Y  if  K  and  F  are  properly  defined.  The  process 
by  which  jC  and  F  are  assembled  and  the  relation  between  Y  and  the 
desired  solution  y(x,t)  will  be  described  here  in  this  section. 

The  first  step  is  to  write  down  the  expressions  in  the  vari¬ 
ational  statement  in  terms  of  the  element  variables.  A  grid  scheme  of 
elements  is  shown  in  Fig.  2.  In  this  figure,  the  nondimens ional  length 
of  the  gun  tube  is  divided  into  K  equal  segments  and  the  time  range  of 
interest  into  L  equal  segments.  The  result  is  then  a  set  of  K  x  L 
rectangular  elements.  In  the  equations  that  follow  n,  the  sub-  or 
super-scripts  m,n  denote  the  association  with  the  mt*1,nt^  segments  or 
the  (m.n)**1  element.  Define  the  relation  between  the  local  coordinates 
(£,ri)  of  the  (m,n) element  and  the  global  coordinates  (x,t)  by 

K  ■  ■  Kx  -  m  +  1 

q  =  r/n)  «  Lt  -  n  +  1  (9) 

0r»  ,  i 

x  =  -  (S>m-1)  ;  t  =  -  (n+n-1) 

K  L  _ _  _  .  . 

A cccTsrl  or  v  •  r 

I  -• 
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Fig.  2  Finite  element  grid  scheme  in  space-  and  time-coordinates. 
Also  shown:  node  point  ordering  in  a  typical  element  (m,n). 


One  then  can  write  from  Eqs.  (S)  and  (6) 

K  L  „3  1  1 

xt  _  r  r  K4  t  t 


«ia  I  I  r  /  I  Y'l  nJy*f''  .mn 
m=l  n=l  L  0  0  lra*ni  (m»n) 

i  1  1 

=  II  f  (p"«  sin  “)/  /  [ (K-m+1)  -  Qy'f  fiy*'  d?dn 


0  0 
,  1  1 


(m,n)  (m,n) 


«3  ‘  -  H  W  / 

514  ■- Ht/ /V;..n)  <w)dEdn 

515  ■  n  “  /  tt-odom  (10.) 


( ) 
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^,*Y  /  /  [(n-l)l+2(m-i)n+n2]y»  Sy*’  6  fE-ndfdn 
L  o  0  C®.n)Oy0n,n)6(m,n)(C  c)d?dn 


*t  rr  p»62Ky2  **  .* 

L7  LI  .]  j  j  £ C*n-1) z-*-2 Cm-l)n+n2ly •  Sy*'  S’  rr 

1  o  0  Jy0n.n)Oy(m,n)6(m,n)CC  C)d?dri 

*'•  '  niSLri/0  0  (5-OdCdri 


fiI9  “  "  II  /  /  yf  .6y*  6  .(C-C)dCdTi 

K  0  0  (m»nJ  (™»n)  (m,n)  ^  ^  ^  1 

dI10  -  -  II  f1  f1  Y,  6y*  i  (i-E)dEdn 

K  0  0  (m*n)  Cm.n)  (m,n)  5,d5dT) 

L  1  1 

8,11  ’  n-1  *  k2lty(l,i»<yU.ii) 

*  Vff.,,)*’*...)  *  Wn/^.n)1'"' 

Ml2  *  J.  "T 

K  L  »  11 

<5 j  _  _  V  r  Y  g  cos  a  r  .  „ 

i  KI  /  /  5y*  d^dn 

m»l  n=l  KL  0  0  (n»n) 

*.t  -  VV  Y2ffn  cos  a  tl  A  .  - 

2  it  “l - K  '  6yr. -A  ,CC-C)dCdn 

L  0  0  (*»n)  (m,n) 

K  ky  1 

6j3  *  I  t  /  Ytm\M*y*tm 

m=l  o  W  (m»L) 


(10b) 


W?) 
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Now,  the  shape  function  vector  a(g,n)  is  introduced 

m.n) 


a  (5,U)x 


(m,n) 


(C,T))  =  a  (C»T))Y* 


(12) 


where  X(m,n)  and  X*(m,n)  are  the  generalized  coordinates  for  the  orig- 
inal adjoint  field  variables.  Thus,  equations  in  the  form  of  Eq.  (7) 
and  (8)  are  obtained  in  the  usual  manner  of  finite  element  formulation 
(see,  for  example.  Ref.  [5)). 


5.  NUMERICAL  RESULTS 

In  actual  computation,  a  simplified  model  of  a  M68-105  mm 
cannon  tube  is  used.  The  source  of  data  can  be  found  in  Reference  [4], 
For  the  cannon  tube,  we  take 

E  =  3  x  107  psi  ;  p  =  0.283  lb/ in. 3 

l  =  210  in.  ;  A  *  28.2  in.2  ;  I  =  123.6  in.4 

and,  for  the  projectile, 

nip  *  20  lb  ;  p  (inbore  pressure)  =  40,000  psi 

$  (projectile  acceleration)  =  6,000,000  in. /sec2 

v  =  25,000  ft/sec  ;  T  =  0.008  sec 

average 

Consequently, 

c  =  (pAJl7(EI))1/2  -  0.104  sec 
A/c  =  2100  in. /sec  ;  y  -  c/T  =  13.00 

In  the  finite  element  discretization,  we  have  chosen  a  grid 
scheme  of  5  x  5  =  25  elements  as  shown  in  Fig.  2.  The  shape  functions 
used  are  cubic  in  x  and  linear  in  t.  As  an  illustrative  example,  the 
support  conditions  are  those  of  a  cantilever.  Thus  k.  and  kj  are 
infinite.  In  actual  computations,  they  are  taken  to  be  very  large 
numbers  compared  with  unity  (e.g.,  108).  The  initial  conditions  are 
such  that  the  tube  has  zero  deflection  and  zero  velocity  at  t  ■  0. 
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TABLE  1.  DEFLECTION  OF  A  GUN  TUBE  [y(x,t)/£]  x  104  AS  THE 
PROJECTILE  MOVES  THROUGH  THE  TUBE  TILL  EJECTION 


0 

0.2 

0.4 

0.6 

0.8 

1.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0151 

0.0192 

0.0182 

0.0184 

0.0281 

0.0000 

0.0413 

0.0786 

0.0735 

0.0783 

0.0986 

0.0000 

0.0724 

0. 1527 

0.1768 

0.1830 

0.2156 

0.0000 

0.0905 

0.2386 

0.3245 

0.3411 

0.3362 

TABLE  2.  SLOPE  OF  A  GUN  TUBE  y'(x,t)  x  104  AS  THE  PROJECTILE 
MOVES  THROUGH  THE  TUBE  TILL  EJECTION 
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Tables  1  and  2  contain  the  results  of  tube  deflection  and 
slope  as  the  projectile  moves  down  the  tube  from  the  breech  end  to 
the  muzzie.  For  example,  at  t  *  0.5t,  the  deflection  represents  tube 
deflection  when  the  projectile  is  at  a  location  of  x  =  0.25Jlfrom  the 
breech  end. 

Figure  3  depicts  the  deflection  curve  of  the  gun  tube  at 
projectile  ejection  and  Figure  4  illustrates  the  muzzle  motion  as  the 
projectile  moves  down  the  tube. 

From  Tables  1  and  2  it  is  observed  that  the  initial  data 
have  been  recovered  for  this  present  unconstrained  variational  formu¬ 
lation.  The  data  presented  here  are  of  the  same  order  of  magnitude 
compared  with  other  analysts  [4].  Further  verification  of  the  compu¬ 
tational  results  and  a  parametric  study  of  the  relative  importance 
of  the  various  forces  involved  in  the  problem  will  follow  in  the 
immediate  future. 
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